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Abstract 

We calculate Df — > <j> transition form factors V, Aq, A\ and A2, and 
study semileptonic decay of Df — > 4>lv based on QCD sum rule method. We 
compare our results of the ratios of V(0)/Ai(0), A2{<$) / A\{Q) , Fl/Tt, and the 
total decay branching ratio of — > (piu with experimental data, and find 
that they are consistent. 

1. Introduction 

Semileptonic decay of charm meson is important for studying strong and weak in- 
teractions. It can be used to test techniques developed for solving perturbative 
and nonperturbative problems in Quantum Chromodynamics (QCD), and to ex- 
tract elements of Cabibbo-Kobayashi-Maskawa (CKM) matrix. Semileptonic decay 
is simpler than hadronic decay of charm meson because leptons do not involve strong 
interaction. The amplitude of semileptonic decay can be decomposed into several 
transition form factors due to Lorentz property of the hadronic matrix element. The 
form factors include all the nonperturbative effects. Several methods can be used to 
treat these problems, such as quark model, QCD sum rule and Lattice QCD, among 
which, QCD sum rule and Lattice QCD are based on the first principle of QCD. 

The method of QCD sum rules jTj has been widely used in hadronic physics since 
its establishment in the late 1970s. For semileptonic decays of charm meson, D + — > 
K e + v e was firstly studied in QCD sum rule method with three-point correlation 
function Several years later, QCD sum rule method was extended to semileptonic 
decays of B meson, B — > D(D*)£i> [3] and B — > nets In these works, form factors 
f+{q 2 ) and fv(q 2 ) are calculated at the point q 2 = 0, where q 2 is the momentum 
transfer squared. For the whole physical region of < q 2 < q^ax, the form factors 

2 

are either assumed to be pole dominance /(0)/(l %— ), or a linear approximation 
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was used. In Refs E], D — > K°e + u e , K° e + u e and D — > 7rei/, pei/ were studied, 
where QCD sum rule method was extended to a very large value of g 2 with a careful 
treatment of non-Landau-type singularities. D s decays to 77 and rf final states were 
studied in [7j. 

In this work, we study — > <^i/ in QCD sum rule method. This decay mode 
has been measured in experiment long time before [HI El EES HI] • Now It is necessary 
to analyze it theoretically. We calculate up to contributions of operators of dimen- 
sion 6 in the operator product expansion (OPE) and keep the mass of s-quark. In 
our result, the large contributions come from unit operator I (result of perturbative 
diagram ) and condensate of operators of dimension 3. Operator of dimension 5 
gives smaller contribution. The contributions of operators of dimension 6 are negli- 
gible. When calculating contribution of perturbative diagram and gluon condensate 
(operator of dimension 4), Cutcosky's rule has been used. Therefore subtraction of 
continuum contribution is conveniently performed not only for perturbative diagram 
but also for contribution of gluon condensate. After some long steps of calculation, 
we find the contributions of diagrams for the gluon condensate cancel each other, 
so there is no gluon-condensate contribution in D£ — > <fi transition. This is our new 
finding. 

There are two independent Borel parameters Mf and Mf in manipulating three- 
point correlation functions. In general, to simplify the numerical analysis, a fixed 
ratio of M x 2 /M| was taken in recent references. In this work, we make the numer- 
ical analysis in the whole region of independent Mf and Mf to select the stability 
"window" , so it is different from taking a fixed ratio of these two Borel parameters. 

We calculate D s — > transition form factors V, A , A\, A 2 and the branching 
ratio of D+ -> <f>£u. Our result of the ratios of V(0)/Ai(0), A 2 (0)/A 1 (0), V L /Y T 
{Tl and IV denote decay width of to <fi meson in longitudinal and transverse 
polarization, respectively), and the total branching fraction are in agreement with 
experimental data. 

Recently, just before this work is finished, we find that — ► <f)£i> was also 
calculated in Ref . [12] . However, their analysis is very different from ours. First, by 
carefully choosing the requirement that the double Borel parameters Mf and Mf 
should not be too large for keeping the continuum contribution small, and at the 
same time, Mf and Mf should not be too small for keeping the truncated OPE series 
effective, i.e., keeping the contributions of higher dimension operators small, we get 
very different stability "window" for the Borel parameters. Second, our results of the 
transition form factors are different from theirs. Especially for A 2 , they got negative 
value, however, we get positive. Using their values of the form factors, although one 
can get the total branching ratio of —>■ <f)£i> to be compatible with experimental 
result, the ratio of Tl/Tt will be too large. But in our case, Tl/Tt = 0.99 ± 0.43 
which is consistent with world average 0.72 ± 0.18. 

The paper is organized as follows. In section 2, we briefly introduce the QCD 
sum rule method used in this work. Section 3 is the calculation. Section 4 is the 
numerical analysis and discussion. Section 5 is devoted to the summary. 
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2. The method 



To calculate the transition form factors of semileptonic D s meson decays, the stan- 
dard procedure in QCD sum rule method is to consider the three-point correlation 
function defined as 

IV = ? / ^Ae^--^^(0|T{^(x)^(0)j 5 D ( 2 /)}|0), (1) 

with the currents having the same quantum numbers as the relevant mesonic states 
under consideration, which are defined by: 1) the current of D s channel, j^{y) = 
c{y)i^s{y)] 2) the current of weak transition: j M (0) = S7 M (1 — 7s)c; 3) the current 
of (p channel: j$(x) = s(x) r Y u s(x). On one hand, inserting a complete set of interme- 
diate hadronic states into the correlation function, and using the double dispersion 
relation, one can express the correlation function in terms of a set of hadronic states, 



,, / , , p(si,s 2 ,q 2 ) 
11 ^ = / dsids 2 - f 577 o\, [2) 

1 (si - pi ){S2 - m 



with 

p(su 82, q 2 ) = Y,(Q\jt\X)(X\j,\Y)(Y\j°\0)5( Sl - m 2 Y )5(s 2 - m 2 x )9(p x )9(p Y ), 

XY 

where X and Y denote the complete set of hadronic states of and D s channels, 
respectively, px and py are the four-momentum of X and Y states, s\ = p Y , 
s 2 = p 2 x, and q = p\ — p 2 . Integrate over s\ and s 2 in Eq. (j2j), we can obtain 

n ^ (0\ji\X}(X\j,\Y}(Y\£\0} 

liuy = > -, — o 777 — o ^ 1" continuum states. Id) 

xy ( m Y ~ Pi)( m x ~ Pi) 

Separate the ground states of D s and channels apparently, the above equation 
becomes 

(oi^0)(0ij^mij 5 D io) , , . , , .. m 

llui/ = 7 — 5 ^77 — 5 q7 h higher resonances and continuum states. (4) 

(m 2 Ds -pf){mj-p 2 2 ) 

The weak transition matrix element D s — > <p can be decomposed as 

rriD s + 
£* • q 

~^ £ l -^-%)( m D a + m^Atiq 2 ) (5) 

+i[{Pi + P2)t, 2 — -g^£ ■ g- 



g 2 m Ds + 



-« 2 — <7^o(g 

g 
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where q = p\—p 2 - The vacuum-to-meson transition amplitudes can be parameterized 
through defining the corresponding decay constants, 

{0\sj„s\4) = m^e^, 

(0\si l5 c\D s ) = (6) 
m c + m s 

Finally the correlation function can be expressed in terms of meson decay constants 
and D s — ► <f> transition matrix element, 

n _ rntheW (</>(4 A) , p 2 ) | Jjj, | D s (Pi ) ) fpm 2 Ds 



(m 2 Ds - p?)(m} - p\){m c + m s ) 
+higher resonances and continuum states. (7) 

On the other hand, the correlation function of Eq. can be evaluated at neg- 
ative values of p\ and p\ by the operator-product expansion in QCD, in which the 
time-ordered current operators in Eq. (JIJ is expanded in terms of a series of local 
operators with increasing dimensions, 

z 2 J ^ye^^T{j*{x)3Mi?{y)} 
= C^J + C^M + C A , u G a a(3 G a ^ + C 5 ^a af3 T a G aa ^ 

+ c 6F *rwr'f + -, (8) 

where C^'s are Wilson coefficients, / is the unit operator, is the local Fermion 
field operator of light quarks, G a a ^ is gluon strength tensor, T and T' are the matrices 
appearing in the procedure of calculating the Wilson coefficients. Sandwich the left 
and right hand sides of Eq. (jHJ) between two vacuum states, we get the correlation 
function in terms of Wilson coefficients and condensates of local operators, 

IV = »^/d 4 a;^e^^-»(0|r{;*( a :)j M (0)jf(y)}|0> 

= C Q ^J + C 3ab/ (0|M|0) + C ilMU (0\G a a(3 G aa ^\0) + C^(0\^a aP T a G aa ^\0) 

+ c 6Att/ (o|*rwr'*|o) + ..., (9) 

For later convenience, we shall reexpress the above equation. In general, it can be 
expressed in terms of six independent Lorentz structures 

Each fi includes perturbative and condensate contributions 

f> = fT + fl® + fl A) +fP + fi* i + (ii) 

where ^ , ■ • are contributions of condensates of dimension 3, 4, 5, 6, ■ • • 
in Eq. Q. In next section we can see that perturbative contribution and gluon 
condensate contribution can be finally written in the from of dispersion integration, 

cpert 



fr = I dSldS2 \s,-Pi)(s,-pir 

,(4) f, , pf\s ll s 2) q 2 ) 

u ~ r^^-pi^-piy 
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We approximate the contribution of higher resonances and continuum states as 
integrations over some thresholds and s® in the above equations. Then equate 
the two representations of the correlation function in Eq. (JJJ) and (|1U|). we can 
get an equation for the form factors. To improve such equation, we make Borel 
transformation over p\ and p\ in both sides, which can further suppress higher 
resonance contribution. The definition of Borel transformation to any function f(p 2 ) 
is 

p 2 — > — OO 
-p 2 /n = M 2 

Some examples of Borel transformation is given in the following, 

B, 1 = 1 1 e -*/M* 

I P 2 > M2 ( s -p 2 ) k (k-l)\(M 2 ) k 
£| p 2 >M2 (p 2 ) fe = 0, for any k > 0. 

Equating the two representations of the correlation function, subtracting the 
higher resonances and continuum contribution, and performing Borel transformation 
in both variables p\ and p\ we finally obtain the sum rules for the form factors, 

2m^f Ds m z Ds 

A^q 2 ) = - K + ( ms) -^s/^e m V^MlMlBh, 

m rt> f^f D m 2 Ds (m Ds + m^) 

A 2 (q 2 ) = K + m)(m D + ^/M^^/Mj^ \ ^ + ^ (12) 
m<f>f<j)jD a m Ds 2 



(m c + m s ) a /M a _ m3 /M 2 2 - 6 , x m D s - m 
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A)(g 2 ) = _ ^ I e ^J^e m V^M 2 Ml[B{h + / 3 
+fi(/i-/ 3 )y +/b], 

where S/j denotes Borel transforming /j in both variables and p|, M\ and M2 are 
Borel parameters. Because we have subtracted the higher resonance and continuum 
contribution, now the dispersion integration for perturbative and gluon condensate 
contribution should be performed under the threshold, 

ir = r°'d Sl r l i S2 «..«■) 



{si-p\){s 2 -piy 



(si-pf)(s 2 -p|) 



In the next section, we will explain the technique of calculating the Wilson 
coefficients and given the resulted form of the sum rules for the form factors. 
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3. The Calculation of the Wilson Coefficients 



In this work, we first calculate the Wilson coefficients in the operator-product ex- 
pansion jTSJ , then extract the relevant terms fi 's for the sum rules of the form factors 
in Eq. (|12jl . We will not present the result of each Wilson coefficient here because 
their forms are very tedious. We only give the results of the form factors according 
to the contribution of each condensate. 

3.1 The Calculation of the Perturbative part 

The diagram for the perturbative contribution is depicted in Fig.l. The leading 
order in a s expansion is considered here. This contribution amounts to Wilson 
coefficient Cq in OPE representation of the correlation function in Eq. Q. We can 
write down this amplitude (see Fig.l), 




Figure 1: Diagram for perturbative contribution. 



Co 



d 4 k 



l)Tr 



(2?r) 4 

(k + pi) 2 — m 2 + is 



l ^T~2 TV~^ u Ti — I \a 2~< ~7m(1 _ 7b) 

k — rrig + is \k + P2) — vn A s + %e 



(13) 



The above integration can be performed according to Cutkosky's rule ^4j. That is, 
to write the integration of Eq. ()13|) in the form of dispersion integration, 



ds\ds2- 



si -n){s 2 -pi) 



(14) 



The spectral density p(si, s 2 ,q 2 ) can be directly calculated by substituting the de- 
nominators of the quark propagators for S functions, i.e., putting all the quark lines 
on-mass-shell, 

— -> -2iii5{k 2 - m 2 s ), etc., (15) 

k 2 -m 2 + ie 



then the spectral density can be calculated from, 



p(si, s 2 ,q 2 



'-2m) 



d A k 



(27ri) 2V ' J (2tt) 4 

x(fc+ t>i + m c )}5{k 2 - m 2 )5[{k + p x ) 2 - m\\ 
x5[(k + p 2 ) 2 -m 2 2 ] 



Tr[j 5 (fc + m s )j u ( i fc+ p 2 + ra^T/X 1 - 7s) 

(16) 



Pi -> si 
p\ -> s 2 



To perform the above integration, some basic formulas are needed. Part of them 
have been given in Ref.p3] without the quark mass, here we calculate them with 
the quark masses included, 



I = J d A kb{k 2 - m 2 )S[(k + Pl ) 2 - mj)5[(k + p 2 ) 2 - m 
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2VA' 



d^kk^k 2 - m 2 )5[{k + Pl ) 2 - m{]8 [(k + p 2 ) 2 - m 2 ] 



= aipxp + b x p 2fl 



-5^72 [ s 2(-si + s 2 - q 2 ) + (si + s 2 - q 2 ){m 2 - mj) 

— 2s 2 (m 2 — ml)}, 

-^72 [si(s 2 + si- q 2 ) + (si + s 2 - q 2 )(m 2 - m\) 

— 2si(m 2 — ml)}, 



(17) 



;is) 



1^ = I d 4 kk^k u 5{k 2 - m 2 )5[{k + Pl f - m\}5 [(k + p 2 ) 2 - ml] 
= a 2 p ltl pi u + b 2 p 2iiV + c 2 {p XjX p 2v + p\ v p 2p ) + d 2 g^ v , 



a 2 
b 2 
c 2 

do 



s\ — m\ + m 2 , 



D 2 = s 2 - ml + m 2 , 



^m 2 s 2 + {[3s 2 D iai - ( Sl + s 2 - g 2 )£> 2 &i + s 2 D 2 b x ] 
■jj I ^m 2 s 1 + \[siD x ai - {s x + s 2 - q 2 )D x bi + 3si-D 2 &i] 



•^m 2 §(si + s 2 - g 2 ) 



+ s 2 - g 2 )Aai - + |(si + s 2 - g 2 )/J 2 6i] 

_^ + l [jDiai + jDA ] 



(19) 



where A(s 1; s 2 , q 2 ) = (si + s 2 — g 2 ) 2 — 4siS 2 , and in Eqs. (jT7jl ~ substitutions 
p 2 — ► Si and p 2 — > s 2 have been indicated. 

3.2 The Contribution of bi-quark Operators ^(0)^(x) 

The diagrams for the contributions of ^f(x)^f(y) and 1 l , (0)\E'(x) are shown in Fig. 
2. The contribution of Fig. 2(b) is zero after double Borel transformation in both 
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variables p\ and p\ because only one variable appears in the denominator 1 / (p 2 , — m 2 s ) . 
So we will not consider Fig. 2(b) in the following. The contribution of Fig. 2(a) to 
the correlation function is 




Figure 2: Diagrams for the contributions of non-local operators *$>(x)*$>(y) and ^(0)^(x). 



n£ = i 2 J dVye^-^-y(o\y( x ) lu tS s F (xh»(i - K)iS c F (-y)ns*(y)\o), (20) 

where iS F (x) and iSp(-y) are the propagators of s and c quarks, respectively. Move 
the quark field operators *&(x) and ^f(y) together, we get 

ti 2 ; u = i 2 j ^V P2 ^ ipi - y ^ 

(21) 

where a and f3 are Dirac spinor indices. The matrix element (0\ty p(x)ty a (y)\0) can 
be dealt with in the fixed-point gauge JB]. We expand it up to the order of x 3 and 
y 3 using the technique explained in ITS] ITTj. 



(0|f-(rr)^(y)|0> = 5 ab 



? T77^ \ — / 1 

"3! 96 (x " y)2( ^ + ^ TG(7 *> (tq^^-^V 

+ •••], (22) 

a and 6 in the above are the color indices, m is the quark mass, and the ellipsis 
stands for terms of higher orders in x and y expansion. From Eq. fl22(l we know that 
Fig. 2 (a) contributes to the coefficients of quark condensate (^f^f), mixed quark-gluon 
condensate g(^/TGa^/) and the four-quark condensate fifty) 2 . Substitute Eq. (J22j) 
into ()21|) and integrate over the coordinates x and y, we can obtain explicitly the 
coefficients of these condensates contributed by Fig. 2 (a). 
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3.3 Contributions of bi-Gluon Operator G a ^G a ^ 

The diagrams for the contribution of bi-gluon operator are depicted in Fig. 3. They 
are calculated in the fixed-point gauge, in which the gauge fixing condition is taken 
to be x^A^x) = [16 . Then the external gauge field can be expressed directly in 
terms of the color field strength tensor 



daax p G a nn {ax) 



(23) 



which is expanded to the first order to be 

1 



(24) 



In the following calculation, it is convenient to transform A a (x) to the momentum 

space, 

' ^2 d xAfu ^ a , n , , (25) 



Al{k) = --{2irYi r 5\k)G a p ^) + 



dk n 






(b) 



(c) 






(d) 



(e) 



(f) 



Figure 3: Diagrams for contributions of bi-gluon operator. 



Then the amplitude can be written down in the momentum space by following 
the standard Feynman rule. Again, as what we did in the previous subsections, 
we move the gluon strength tensor operator together: G a aa G h ^ p . Then using the 
following decomposition to obtain the bi-gluon condensate, 



(Q\G a arj G b Pp \o) = ^(GG)5 ab (g aP g ap - g ap g a p): 



(26) 
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in which (GG) is the abbreviation of (0|G£ 1/ G WH/ |0). 

In the evaluation of the diagrams of Fig. 3 some types of loop integrals encoun- 
tered are treated at first by derivatives with respect to the quark masses, then 
transform them to dispersion integrals by using Cutkosky's rule and with help of /, 
1^ and functions given previously. For instance, 

d A k- M '' 



(k 2 - m 2 ) [(k + p 2 ) 2 - m 2 } 2 [{k + Pl ) 2 - m\\ 2 
d r k,,ku 



and 



dm 2 , dm\ J {k 2 — m 2 ) [{k + p 2 ) 2 — m|] [(k + P\) 2 — ml] 
d d f I 

= ~ 27ii Tr^Tr^i / ds i ds 2- ( 2T7 2v ( 27 ) 

(77712 V m l J ( s l — Pl) ( s 2 — P2) 



d 4 k k^kyk ■ P2 



{k 2 - m 2 ) 2 [(k + p 2 ) 2 - ml] 2 [(k + pi) 2 - m\\ 2 
d d k ll k v k-p 2 



dm\ dm 2 J (k 2 — m 2 ) 2 [(k + p 2 ) 2 — m 2 .] [(k + Pi) 2 — m 2 } 

_ 2ni A ? ^_ r dsidS2 zMf? - m l + m2 ) V , 28 ^ 

dm 2 dm\dm\j (s\ — p\){s2 — pl) 

where the term — |(s 2 — 1^2 + 'm 2 ) comes from the S functions 8{k 2 — m 2 )5[(k+p2) 2 — 
m 2 ,] with the substitution p\ — > $2 when using the Cutkosky's rule. 

After some long steps of calculation, we finally find that the contributions of the 
diagrams (a) to (f) in Fig. 3 cancel each other. Therefore there are no contributions 
of gluon condensate in D s — > $ transition. 



3.4 Contributions of Quark-Gluon mixing and Four-Quark 
Operators: ^{x)^{y)G^ v and (W) 2 

The diagrams for quark-gluon mixing and four-quark contributions are depicted 
in Fig. 4 and Fig. 5, respectively. The techniques are similar to that explained in 
previous subsections. We only give some different points here. 

The vacuum average of non-local quark-gluon mixing operator ^ {x)^ {y)G a is 
calculated to be 



<0|^ a (s)ttJ(y)GJUO) 
-L(y (7 TG*)(a» l/ ) Pa T« + 



(J {^) 2 (g p ^-g p ^)(x + yY 



96 x 9 

T£, (29) 

/3a 



where (tyaTG^) and (##) 2 are the abbreviations of (0|*<7 Mi ,T a G a ^|0) and (0|^|0) ; 
respectively, g is the strong coupling. 



10 



(a) (b) (c) (d) 

Figure 4: Diagrams for mixed quark-gluon operators. 






(b) 



\ / 



(c) 






(d) 



(e) 



(f) 



Figure 5: diagrams for four-quark contributions. 



Because we calculate up to the condensate of dimension-six operators, the ex- 
ternal gluon field A a Jx) in Fig. 4 should be expanded up to the second term, which 
will contribute a dimension-six operator, 



A^x) 



daax p G a pp {ax) 



^G a pp (0) + ±x«x p D a G a pp (0) + 



(30) 



where D a is the covariant derivative in the adjoint representation, (D a ) mn = d a 5 r 
gf amn A a a . Then another vacuum matrix element needed is jTH] 
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■(m) 2 (g^ a -g^7 P )f3 a T« . 



(31) 



. P±^£\J ap \„ ! 33 X 2 4 N " " ' y ^"' P ' a IP/Put^ Jl 

We calculate these diagrams and find that the contributions of Fig. 4(c), (d) and 
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Fig.5(c), (d) vanish after double Borel transformation in two variables p\ and p|, 
because only one variable appearing in the denominator, for instance, g 2( p /„ m ^) ■ The 
Borel transformation in pf will kill such terms. 

Following the above method, after some tedious algebraic derivation with the 
software MATHEMATICA, we obtain the coefficients f , fi + / 3 , fi — f 3 and f 5 
needed in Eq. (fl2"|). They are listed in Appendix. 

4. Numerical Analysis and Discussion 

In the numerical analysis the standard values of the condensates at the renormal- 
ization point \i = lGeV are taken [T| |2"U]. 

(qq) = -(0.24 ± O.OlGeV) 3 , (ss) = m 2 Q (qq), 
g^aTm) = m§(*tf), a s (M) 2 = 6.0 x 10- 5 GeV 6 , (32) 
ml = 0.8 ± 0.2GeV 2 . 

The quark masses are fixed to be m s = 140MeV, m c = 1.3GeV [21] . and the decay 
constant of <ft meson is extracted from experimental data = 0.228 [22]- For the 
decay constant of D s meson we take fr> s = 0.214 ± 0.038GeV [2*Tj . 

The Borel parameters M\ and M 2 are not physical parameters. The physical 
result should not depend on them if the operator product expansion can be calcu- 
lated up to infinite order. However, OPE has to be truncated to some finite orders 
in practice. Therefore, Borel parameters have to be selected in some "windows" to 
get the best stability of the physical results. The requirement to select the stable 
"windows" is: the Borel parameters can not be too large, or, contributions of higher 
resonance and continuum states can not be effectively suppressed; at the same time, 
they should not be too small, or, the truncated OPE would fail because the series in 
OPE generally depend on Borel parameters in the denominator 1/M. We find the 
optimal stability with the requirements shown in Tableland the thresholds s°, s° in 
the ranges = 5.8 — 6.2GeV 2 , s° = 1.9 — 2.1GeV 2 . The regions of Borel parameters 
which satisfies the requirements of Table ^ are shown in FigJB] in two-dimensional 
diagram of M\ and Mf . We find good stability of the form factors within these 
regions. 

Because it is not easy to show the contribution of each term of OPE in two- 
dimensional regions of M\ and Mf, we show the contributions of perturbative and 
condensate terms in Table[2]at a representative point (M 2 , Mf) in the stable region 
of Mf and Mf . In general the higher the dimension of the operators, the smaller 
the relevant contributions of the condensates. The main contributions to V(0), 
Ai(0) and A 2 (0) are from perturbative and quark condensate term. For A (0), 
the largest two contributions are from perturbative term and mixed quark-gluon 
condensate. Contributions of four-quark condensate are less than a few percent, 
therefore contributions of operator of dimension 6 are negligible. 
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Jiff (GeV 2 ) 



Mf (GeV 2 ) 



! V 



1.5 2 2.5 3 3.5 4 4.5 5 



Mf (GeV 2 ) 



M 2 (GeV 2 ) 




1 1.5 2 2.5 3 3.5 



00 



M 2 (GeV 2 



Mf (GeV 




M 2 (GeV 2 



(c) 



(d) 



Figure 6: Selected regions of Mf and Mf: (a) for V; (b) for A ; (c) for A\; (d) for ^4 2 - 
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Table 1: Requirements to select Borel Parameters M 2 and M| for each form factors 
V(0), A (0), A x (0) and A 2 (0) 



Form Factors 


contribution 


continuum of 


continuum of 




of condensate 


D s channel 


cf) channel 


no) 


<49% 


<5% 


< 26% 


4)(o) 


< 29% 


< 16% 


<31% 


Mo) 


<49% 


< 18% 


< 22% 


A 2 (0) 


< 11% 


< 27% 


<5% 



Table 2: Contributions of perturbative and condensate terms in the operator- 
product expansion to the form factors V(0), A (0), Ai(0) and A 2 (0), at a repre- 
sentative point (Mf, Mf) in the stable region of M 2 and Mf. fP ert ; perturbative; 
f^: quark condensate; f^: gluon condensate; f^: mixed quark-gluon condensate; 
f^: four-quark condensate. 



Form Factors 


total 


jpert 


;(3) 


;(4) 


;( 5 ) 


/(«) 


{Ml Ml) GeV 2 


V(0) 


1.20 


0.63 


0.66 





-0.10 


0.01 


(2.2,1.4) 


4,(0) 


0.43 


0.28 


-0.10 





0.23 


0.02 


(1.7,1.5) 


Ai{0) 


0.53 


0.28 


0.20 





0.04 


0.01 


(2.0,1.2) 


A 2 (0) 


0.57 


0.22 


0.44 





-0.09 


0.00 


(3.6,1.5) 



The final results for the form factors at q 2 = are 

V(0) = 1.21 ±0.33, Aq{0) = 0.42 ±0.12, 
Ax(0) = 0.55 ±0.15, A 2 (0) = 0.59 ±0.17, (33) 
V(0) A 2 (0) 

rv =i^j = 2 - 20±0 - 85 - r ^^d = 1 - 07±0 - 43 - 

We compare our results for the ratios of form factors with experimental data in 
Table El It shows that the results are consistent with experimental data. 

The physical region for q 2 in D s — > <fitv decay extends from to {rri£, a — m^) 2 ~ 
0.9 GeV 2 . In the range q 2 < 0.4GeV 2 , there is no non-Landau-type singularity [S] 
with the thresholds and s 2 chosen in this paper. The q 2 dependence of the form 
factors is shown in Fig|7|in the range -0.4GeV 2 < q 2 < 0.4GeV 2 . Within this range, 
the behavior of V(q 2 ) and A Q (q 2 ) is well compatible with the pole-model, 



l-q 2 /m^ ole ' 

While the q 2 dependence of and Ai(q 2 ) and A 2 (q 2 ) is very weak. 
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Table 3: Comparison of our results of ry and r 2 with experimental data: E791 is 
from Ref.dU, CLEO from Ref.pU, E687 from Ref.p and E653 from Ref.0. 





ry 




E791 


2.27 ±0.35 ±0.22 


1.57 ±0.25 ±0.19 


CLEO 


0.9 ±0.6 ±0.3 


1.4 ±0.5 ±0.3 


E687 


1.8 ±0.9 ±0.2 


1.1 ±0.8 ±0.1 


E653 


2.3^ ±0.4 


2.11&jj ± 0.2 


Average 


1.92 ±0.32 


1.60 ±0.24 


our result 


2.20 ±0.85 


1.07 ±0.43 




. 8 
. 6 



0.4 

-0.4 -0.2 6 0?2 0.4 

q 2 GeV 2 



Figure 7: q 2 dependence of the form factors from QCD sum rule. The solid curve is for 
V(q 2 ), the short dashed curve for Ao(q 2 ), the long dashed curve for Ai(q 2 ), and the dotted 
one is for A2{q 2 ). 



We fit V(q 2 ) and A (q 2 ) by the pole model in the range — 0.4GeV 2 < q 2 < 
0.4GeV 2 , and extrapolate the fitted result to the whole physical region. The fitted 
pole masses are, 



v 

m pole 



2.08 ±0.13 GeV, 
1.9 ±0.2 GeV. 



, pole - x.^u.^v. (34) 

The form factors calculated in QCD sum rule in this paper are used to calculate 
the differential and total decay rate of D s — > <piv decay. There are three polarization 
states for meson: one longitudinal state, two transverse polarization states (right- 
handed and left-handed). The differential decay rate to longitudinally polarized <p 
meson is 



~dq~ 2 



G 2 F \V CS \ 2 
1927r 3 m^ 



X(m 2 D .m 2 ,q 2 



2m^ 



(m 2 Dg - mj - q 2 ){m Ds + m^A^q 2 ) 



X(m 2 D ,ml,q 2 



-A 2 (q 2 



m Ds + 

where Gp is Fermi constant, V cs is CKM matrix element for c 

•Li. + m 



(35) 

s transition, and 



\ f 2 2 2\ / 2 i 2 



q 2 ) 2 - Am 2 D ml. 
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The differential decay rate to transverse state is 

V{q 2 ) 



dT$ 
ltd 2 



G F | V cs | 2 

1927T 3 



\ / 2 2 2\ 



in: 



=F 



(mc s + m0)Ai(g 2 ) 



™d 3 + ™4> J\(m? Da ,ml,q2) 



(36) 



where the symbol " + " and " — " denote right and left-handed states, respectively. 
Finally, the combined transverse and total differential decay rates are 



dT 7 
dq 2 



dT 



d 



dq 2 dq 2 



(r L + r T ). 



(37) 




Figure 8: Differential decay widths of Df — ► <piv as a function of momentum transfer 
squared q 2 in unit of 10~ 14 GeV _1 . 

The differential decay widths as a function of momentum transfer squared q 2 are 
shown in FigJHJ Integrate them over q 2 in the whole physical region from q 2 = to 
( m D a — m ti>) 2 , we get the integrated decay widths 

T T = (1.39 ± 0.75) x lCT 15 GeV, = (1.05 ± 0.22) x 10- 14 GeV, 

T L = (1.18 ± 0.43) x 10- 14 GeV, T T = (1.19 ± 0.29) x lQ- 14 GeV, (38) 
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and the ratio of Tl/Tt is 

T L /T T = 0.99 ± 0.43 , (39) 

which is consistent with the averaged experimental data {V l/TtY xp = 0.72 ± 0.18 
|2"2~] . The detailed comparison of this ratio with experimental data is shown in Table 

m 

Table 4: Comparison of our results of Tl/Tt with experimental data: CLEO from 
Ref.pU, E687 from Ref.0 and E653 from Ref.0. 





Tl/Tt 


CLEO 


1.0 ±0.3 ±0.2 


E687 


1.0 ±0.5 ±0.1 


E653 


0.54 ±0.21 ±0.10 


Average 


0.72 ±0.18 


our result 


0.99 ±0.43 



We use the total decay width of D s meson Tn s = 1.34 x 10 12 [22] to obtain the 
branching ratio of Df — ► <fi£i/, our result is 

Br(D+ -> <f£v) = (1.8 ± 0.5)%, (40) 

which is in good agreement with experimental data Br(D+ —>■ <ft£i>) exp = (2.0 ± 
0.5)%. 

5. Summary 

We calculate the transition form factors for D s — > <p transition in the region q 2 < 
0.4GeV 2 in QCD sum rule, where no non-Landau-type singularity occurs. Then fit 
the result from QCD sum rule in this region of momentum transfer, and extrapolate 
it to the whole physical region in the decay — > (f>£v. We treat the two Borel 
parameters M\ and Mf as independent parameters, and select the allowed region 
for M\ and Mf by requiring that the higher resonance and continuum contributions 
in D s and channels are not large, at the same time requiring that the condensate 
of higher dimension operators do not contribute too much. We find good stability 
for the transition form factors V, A Q , Ax and A 2 in the relevant two-dimensional 
regions of M\ and Mf . We obtain the results of the transition form factors V, A Q , 
A\ and A 2 in these regions of M\ and Mf . Our result of the ratios of these form 
factors ry and r 2 are well consistent with experimental data. 

We studied the process — > (f)£v with the form factors calculated from QCD 
sum rule. For the transverse polarization state of the final (f) meson, the rate of 
D s decaying to right-hand state is almost an order smaller than decaying to left- 
hand state. The ratio of Tl/Tt and the branching ratio of — > <f)£v are in good 
agreement with the experimental data within the error bars of both the present 
experimental data and theoretical calculation. 
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Appendix 

Borel transformed Coefficients of perturbative and nonperturbative contributions 
to the transition form factors in Eq. ()12j) are given here. The contributions of 
condensate of dimension- six operator (^/^) 2 are numerically negligible, whereas their 
expressions are more tedious, therefore we do not present all of them here. 



1) Results for / : 



Bfo = Bf% ert + Bf ( * ] + Bf™ + Bf™ 



with 

s o s o 3 e -si/Af 1 2 -s 2 /Af 2 2 

Bfn" = [a 2 2 ds2 f s l dsi ^ — „ „ „ ,„ [2s 2 ^ — 2som s m 2 . — s 2 (2to 2 

JO Jirnl ij s l AMfM^ir 2 X 3 / 2 s (Al) 

+Si — ,s 2 + q 2 )m c + m s (— s| + 2ni 2 s s 2 + sis 2 + q 2 s 2 + A)], 

where A = (si + s 2 — q 2 ) 2 — 4sis 2 . The lower integration limit is determined by the condition 
that all internal quarks are on their mass shell |19j . 

L m c 2 

m z c — q A 

~m 2 /M 2 -m 2 /AI 2 

Bf ° = ~ ' GMfMl [^ M l M 2-ml + q 2 WM + (A{M 2 + M 2 ) 



+q 2 )m 4 M 2 - {Ml + M 2 )m 6 s )M? + M 2 (M 2 + M 2 )m 2 c m 2 s {3M 2 
-m 2 s )M 2 + Mim c m s (-3M?Mi + {Ml + M 2 + q 2 )m 2 s M^ 
-{Ml + Mf )mf)M? - M${M 2 + Mi)m 3 c m^\ x (ss) , 

-m 2 /M 2 -m 2 /M 2 

^° = 12M 8 M 8 t(3M| - MlM 4 + {Ml + M 2 )m\ 

-(3M 2 4 + 5M?M 2 2 )m 2 s )M? - M%{Ml + M^)m 2 c {?,M 2 - m 2 )Ml 
+Mlm c m, s {2M 2 1 - 2M 1 2 M| + {Ml + M 2 )m 2 s )Ml + M${Ml + M%)mlm s 
+q 2 {M?M 2 2 {3M 2 2 - m 2 ) - MlM%m c m s )] x g{saTGs) , 



{A2) 



{A3) 
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81MfMl(m 2 c - q 2 )m 3 

+{Ml + M 2 )m 5 c m 3 s M 4 + M 2 (M 2 + M 2 )m 4 c m 4 M 2 
+M 2 q 4 m 3 (m s M 2 + M 2 m c )M 2 + M 2 m 3 c m 3 s (-13M$ + 2M 2 M 2 
+{Ml + Ml)m 2 s )Ml + M 4 to 2 (-54(-1 + e m i/ M 2)M 2 M$ 
+54M 2 m 2 M 4 - (M 2 + 10M 2 )m 4 M 2 + (M 2 + M 2 )m 6 s ) 
+g 2 (54(-l + e m '/ M 2)MlMf - 54M 4 m 2 M 1 6 - (M 2 + M 2 )m 6 s M 4 
-M 2 (M 2 + M 2 )m c m\M 2 + M 2 (M 4 + 10M 2 M 2 - (2M 2 
+M 2 )m 2 c )m 4 M 2 - M 2 m c {2M 4 - 13M 2 M 2 + (2M 2 
+M 2 )m 2 c )m 3 )}xg 2 (ss) 2 . 

2) Results for /i + / 3 : 

B(h + h) = Bf + ert + B/f + Bff + Bff , 

with 

B.fT rt = JZi ds 2 Q 4M 2 M|7r2A5/2 {-6« a (-«i + ^ + 9>« + 6s 2 (- Sl 
+s 2 + q 2 )m s m 4 + 2s 2 (-4s 2 + 8s 2 s 1 - As 2 + 2q 4 - 6(si - s 2 )m 2 
+A + 2q 2 (3ml + si + s 2 ))ml + 2s 2 m s [Asl - 8s 2Sl + 4s 2 . - 2q 4 
+6(si - s 2 )m 2 - 3A - 2g 2 (3m 2 + Sl + s 2 )]m 2 c + [6(si 
-s 2 )s 2 m 4 + 2(4s| - 8sis?, + 4s 2 s 2 - 2Xs 2 + si\)m 2 s + (si 
-s 2 )s 2 (2s 2 - 4s 2 si + 2s 2 . - A) - 2s 2 g 4 (2m 2 + 2si 
+s 2 ) + g 2 (-6s 2 m 4 - 2(2s2. + 2s x s 2 + A)m 2 + s 2 (2s 2 - 6s 2 si 
+Asi, - A))]m c + m s [2s| - 6sis| + 6s 2 s2. - 3AS 2 . + 6(s 2 
-si)m 4 s 2 - s\s 2 + 3siAs 2 + 2g 4 (2m 2 + 2si + s 2 )s 2 + A 2 
-2(4s^ - 8sisl + 4s 2 s 2 - 4As 2 + siA)m 2 + q 2 {6s 2 m 4 + 2{2s\ 
+2sis 2 + A)m 2 + s 2 (-2sf + 6s 2 si - As\ + 3A))]} , 

-m 2 a /M?-m 2 s /M% 

* f + = QMfMl {[6M 2 M|-3(M 2 -2M 2 + g 2 )™ 2 M 2 4 + (4(M 2 

+M|) + q 2 )m 4 M 2 - (M 2 + M 2 )m 6 s }M 4 + M 2 {M 2 + M 2 )m 2 c m 2 (3M 2 
-m 2 )M 2 - M 2 m c m s (3M 2 M 4 + (M 2 + M 2 - q 2 )m 2 M 2 
+ (M 2 + M 2 )m 4 )M 2 - M 4 (M 2 + M 2 )m 3 c m 3 } x (ss) , 
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p-m t ./m 1 -m s /jvi 2 

Bf + = UMfMl {[(M i + 3M * )M * ~ (M ' + M ° )m ' + ( ™ 2 

+5M 2 M 2 2 )m 2 ]M 1 4 + Mf (M- 2 + M 2 2 )to 2 (3M 2 2 - ™ 2 )M 2 - M 2 2 m c m s 

[2(M 2 + 2M 2 )M 2 + (M 2 + Af 2 )m 2 ]M 2 - M 4 (M 2 + M 2 >3 mj + g 2[ M | ( , 
-3M|)M 4 + Af 4 m c m s M 2 ]} x g(scrTGs) , 

p- m ll M l- m \l M l 2 2 

^ - isgpgra - i>g { - 18 '-' + e "' JU ' )Mf m - M ° + (7M| 

-2Af 2 )m5m3M 4 + M 2 (M 2 + M 2 )™ 4 m 4 M 2 + Mfq 4 m 3 s (Mfm s 
-2M 2 m c )M 2 + M 2 m3m3[-26M 4 + 4M 2 M 2 + {M 1 - 2M 2 )m 2 ] 
Mf + M 4 m 2 [-54(-l + e m '/ M 2)M 2 M 2 6 + 54M 2 m 2 M 4 - 4(M 2 
-2M 2 )m 4 M 2 + (M 2 + M|)mf] + g 2 [54(-l + e m ^)MfAf| 
~54M 2 4 m 2 Mf - (M 2 + M$)m 6 s Mf - M 2 2 (M 2 - 2M 2 2 )m c m^M 1 2 
+M 2 (4(M 4 - 2M 2 M 2 ) - (2M 2 + M 2 )m 2 )m 4 M 2 + M 4 m c 
(-4M 4 + 26M 2 M 2 + (4M 2 ~ 7M 2 »^]} x g 2^2 _ 



3) Results for /i — / 3 : 

- / 3 ) = B£ er * + B/i 3) + S/i 5) + S/i 6) , 

with 

n n _3 e -Sl/M 1 2 -s 2 /M 2 2 

4M 2 M 2 tt 2 A 5 /2 { 6s2 ( Sl + 3s2 9)^c-6s 2 (si 
+3.s 2 - g 2 )m s m 4 + 2s 2 [-4s? - 4s 2 si + 8s 2 , + 2q 4 - 6(si + 3s 2 ) 
to 2 + A + 2q 2 (3m 2 + si - 5s 2 )]™? + 2s 2 m s [4s^ + 4s 2 Si - 8s^ 
-2g 4 + 6(si + 3s 2 )to 2 + A - 2g 2 (3m 2 + si - 5s 2 )]m 2 + [6s 2 (si 
+3s 2 )m 4 + 2(si + 2s 2 )(-4s 2 , + 4sis 2 + A)m 2 + (si - s 2 )s 2 (2s 2 
-2s 2 , - A) + 2s 2 q 4 (-2m 2 - 2si + s 2 ) + q 2 (-6s 2 ™ 4 - 2(-10s 2 , 
+2sis 2 + A)m 2 - s 2 (-2s 2 - 10.s 2 si + As\ + Xj)]m c - m s [2s| 
-2sisf — 2s 2 s 2 — As 2 + 6(si + 3s 2 )m 4 s 2 + 2sfs 2 + siAs 2 
+2< ? 4 (-2m 2 - 2si + s 2 )s 2 - A 2 + 2(-8s| + 4sis?, + 4s^s 2 
+4As 2 + siA)m 2 + <7 2 (-6s 2 m 4 - 2(-10s 2 , + 2sis 2 + A)m 2 
+s 2 (2s 2 + 10s 2 si -4s 2 + A))]} , 
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[3M 2 M 2 - (M 2 + 2M 2 )m 2 ]m 4 + [3M 4 m s Mf - 9M 2 (M 2 + M 2 ) 



m 3 s M% + (Mf + AMlM% + 3M?M$)m^]m 3 c + M 2 [-6M?M$ + (Mf 
+4M$M? + 3M 4 )mf - (6Mf + UM?M$ + 5M 4 M 2 )m 4 
+3(3M?M$ + 2M 4 M 4 )™ 2 ]m 2 + M 2 m s [-15M?M$ + (2M? 
+3M 2 M 2 + M$)m* - (2M| + 13M?M$ + llM 4 M 2 2 )m 4 

+(8M 1 2 M| + llM 4 M 4 )m 2 ]m c + M 4 m 2 [ 3M 2 M 6 _ 4(M 2 + M 2) (A16) 

m 4 M 2 + (M 2 + M|)mf + (2M| + 5M 2 M 4 )m 2 ] + M 2 M 2 Vm 2 

[(m 2 - 3M 2 )M 2 + M 2 m c m s ] + q 2 [(6M 2 M 2 6 - (M 2 + 2M 2 >8 

+ (7M 4 + 5M 2 M 2 )™ 4 - 3(M| + 2M 2 M 4 )m 2 )M 4 + M 2 m 2 m 2 

(3M 2 (2M 2 + Mf) - (2M 2 + 3M 2 )to 2 )M 2 + M|m c m s (— 3Mf M 2 

-(3M 2 + 2M 2 )m 4 + (2M 4 + 9M 2 M 2 )m 2 s )M 2 - M 4 (2M 2 + Mf) 

mj?m;?]} x (ss) , 

^ = 24M 8 M | {^(M? + M 2 )m s m5 _ M 2 (M 2 + M 2 )[3M 2 M 2 

-(M 2 + 2M 2 )m 2 ]m 4 + M 2 m s [(3M 4 + 4M 2 M 2 + M 4 )™ 2 
-2M 2 M 2 (5M 2 + 3M 2 )]m3 + [2M 4 (M 2 + 3M 2 )M 4 + (Mf 
+4Mf M 4 + 3M 4 M 2 )™ 4 - 2(3MfMf + 5M 2 4 M 4 )™ 2 ]m 2 + M\m s 
[AM 2 (M 2 + 3M 2 )M^ + 2(M? + 3M 2 M 2 + M$)m* + (Mf 

-8M 2 M 4 - 13M 4 M 2 )m 2 ]m c + M 4 [4M 2 M 2 6 + (Mj 2 + M|) ( Al8 ) 

m 2 (2M 4 + m 4 ) - (M 4 + 5M 2 M 2 )m 4 ] + g 4 [M 2 (m 2 

-3M 2 )M 4 + M 4 m c m s M 2 ] - g 2 [(2(M 2 + 3M 2 )M 4 + (M 2 + 2M 2 )m 4 

-2(2M 4 + 3M 2 M 2 )m 2 )M? + M 2 m c m s (M% - 10M 2 M 2 

+ (3M 2 + 2M|)m 2 )M 2 + M 2 m 2 ((2M 2 + 3M 2 )m 2 

-3(M 4 + 2M 2 M 2 ))M 2 + M 4 (2M 2 + M 2 )m 3 c m s }} x g(saTGs) . 
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